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SUMMARY 



This note extends the method of N.A.C.A. Report No. 
516 to the case of an arbitrary "body of revolution in 
steady motion at constant angle of attack and constant 
path curvature, viz, circular path. Expressions are given 
for the pressure, for the lateral and longitudinal compo- 
nents of the centrifugal force of the apparent mass, and 
for the yawing moment. These expressions are then applied 
to an ellipsoid of revolution mainly for the purpose of. 
demonstrating the usefulness and simplicity of the_y_ector 
method. 

INTRODUCTION 



The circular motion for airship-like bodies has thus 
far "been calculated only for a prolate ellipsoid of revo- 
lution (reference 1, p. 133 and "reference 2), In this 
paper, however, the circular motion of elongated bodies of 
revolution more nearly resembling airships will be inves- 
tigated. The results will give the effect of rotation on 
the pressure distribution and thus yield some information 
as to the stresses set up in an airship in circular flight 

THE BOUNDARY CONDITIONS 



The body is referred to a rectangular Cartesian frame 
OXYZ attached to it, and its axis of symmetry OX moves 
in a plane X'Y' fixed in space. The "body is assumed to 
rotate uniformly in a counterclockwise manner about O'Z' 
with angular velocity CI and to possess an angle of yaw 
P with respect to the origin 0.. (See fig. 1.) 

When considering elongated "bodies, it is convenient • 
to introduce in any meridian plane XH confocal elliptic 
coordinates (j, , A.. The coordinates x, y, z, of any 



2 
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point P in space are then given in this new orthogonal 
system of coordinates by the following expressions: 



x = 2aX\s. 

y = 2a (A. 2 - 1)^ (1 - \f)^'cos co 



= 2a (X s - 1)2 (1 - p. 2 )i 



sin co 



(1) 



The coordinate surfaces of the so-called "spheroidal coor- 
dinates" \i. , X, (0 are obtained by setting, in turn, \x 
constant 1 (hyperboloids of two sheets), X constant (pro- 
late ellipsoids of revolution), co cbnstant (half planes 
through the axis of symmetry OX); these coordinate sur- 
faces furnish by their intersections the three coordinate 
lines:'- \x variable (ellipses), X variable (hyperbolas), 
and co variable (circles with centers bn the OX axis) • 
The coordinate lines give by means of their tangents, di- 
rected positively in the direction in which the corre- 
sponding coordinate is increasing, a rectangular Cartesian 
frame whose origin is the point in which the three coordi- 
nate surfaces (or coordinate lines) intersect. In order 
that this frame be a right-handed one it is so arranged 
that the positive direction of the tangents to the \x, X, 
(0 coordinate lines are analogous, respectively, to OX, 
OY , and 0 2. "■ 



The position of the moving axes is defined by the po- 
sition of the origin 0, whose coordinates with respect 
to the' fixed axes are R cos fit, & sin fit, 0, and 
by the direction cosines of one i*et of axes with respect 
to the other. The following direction cosine table gives 
these values at any time t: 





X 


y 


z 


X' 


sin (Qt---£) 


cos (fit - p) 


0 


■ y' 1 


-cbs (at - p) 


sin (fit - p) 


0 


z \ 


0 


0 


1 



The 'velocity of any point fixed with^ regard to the 
moving axes- is the resultant of two vectors, one of whdch 
V is the Bame for all points of the system, being inde- 
pendent of the coordinates x, y, z of the point and hav- 
ing the components in the directions O'X', O'Y', O'Z 1 
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equal to - R Q sin O.-t, R cos Q't , 0, respectively, 
and in the directions OX, 0Y , 0Z equal to 



= - E Q cos p 
"Vy = - 1 () sin p 
v z = 0 



C2) 



This part of the motion is therefore a translation and JLs 
made up of an axial motion with velocity - R Q cos £ 
and a transverse motion with velocity - R Q sin j3. The 
other part of the motion has components with respect to 
the instantaneous position of the OXYZ frame, given by 



x = - Qy 



(3) 



and is therefore the vector product of a vector H whose 
components with regard to the moving axes are 0, 0, Cl t 
.and of the radius vector F from the origin 6 to the . 
point P. The vector t lies in the XY plane and is 
perpendicular to "both H and r and has a magnitude equal 
to Qr sin (Qr) , It accordingly represents a motion due 
to a rotation of the "body with angular velocity fi about 
the 0Z axis. Referred to a meridian plane making an an- 
gle CO with the XY plane, the vector If has a compo- 
nent t x = - Qh cos CO in the direction of the X axis, 
a component t^ = fix cos (0 in the direction of the H 
axis, and a component Qx sin co perpendicular to the me- 
ridian plane. This latter component plays no part in fix- 
ing the boundary conditions for a body of revolution. 



The translat ional velocity V of the body gives rise 
to motions in the fluid represented by two velocity poten- 
tials: <Pi, due to the axial velocity V x and cp 3 , due 

to the transverse velocity Vy. This type of motion has 
already been discussed in detail in reference 3. Disre- 
garding for the present the translat ional velocity V, 
the normal component of the velocity at any point P of 
the meridian profile, whose element of length is ds, is 
given by 



V n = Q cos to 



(h 



dh 
ds 



+ x 



4s 

ds 



) 



(4) 
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Furthermor-e , if <p 3 represents the velocity pot^en/tlal of 

the flow due to the rotational part of the motion, the nor- 
mal component of the velocity with respect to the Cartesian 
(y,, X, to) frame attached- to the point P, is 

_ ac P 3 _ 3cp a ^X _ gPg a^u, 
9n 3s^ ds Bs^ ds 

a<p . S q> 

or -— ^ ds\ - ~ — ds,. = Q cos co (hdh + xdx) (.5) 

where dsj^, ds^ are the linear elements corresponding to 
the \, p, coordinate lines, respectively. 

Since x = 2a.\\x, h = 2a (\ 2 - l)^ (l - |X S 



and d^ = 2a (^T^ as X = 2a {^ZT^ 



d\ 



it follows that 



(1 ^ 3) SIT -i) 3\ d ^ 

„£ti °* = ( 2a ) n cos a>(|xd|jL+\d\) U> 

C\ 3 - 1) (1 - M- 3 ) " 

Ihe velocity potential cp 3 satisfies Laplace's equa- 
tion and may, in general, he represented by the following 
expansion: 

2 Z C n ra P n m (ix) Q^U) cos HW+-Z ZD/ P n m (n) Q n m (\)sin mo) 

where P n m (p,) , Q n m (X) denote the associated Legendre func- 
tions of degree n and order m of the first and second 
kind, respectively; and C n m , signify certain con-. 

stants which are to "be determined by the boundary condi- 
tions. From the form of the boundary condition (6) it- is 
clear that 

00 

(?) 



cp 3 = Z c a l P n x (p.) Q^ 1 (X) cos CO 
n= i 
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The functions P^ 1 (p,) , Q,^ 1 (A) are given in terms of the 

corresponding Legendre polynomials "by means of the follow- 
ing expressions: 



(8) 



where P n (jx) , Q n (X.) satisfy Legendre' s differential 

equation. T7ith the use of equations (7) and (8) together 
with Legendre' s differential equation for P n (jj,) and 

Qn (X) » *he "boundary condition ( 6) takes the following 
f .o rm: 



CO 

S 
n=i 



a(\|i) dp dQ n a. , 

-dTrdTrdT-- n(n+1) d^n> 



-) 



(9) 



whore it is assumed that A, = X, (p) describes the me- 
ridian curve of the "body in elliptic coordinates. 



For example, if the "body is a prolate ellipsoid of 
revolution its meridian curve is given simply "by X = \ Q , 
a constant. Equation (9) then "becomes: 

and this expression must he valid for the entire range of 
IX. How, of the complete set of P n ' s only one of them, 

namely, P a = \ ( 3p, 2 - l) , .has its derivative proportional 

to |JL. Hence, n = 2 is the only term appearing in the 
"boundary condition for the ellipsoid. Therefore (reference 
1, p. 133) , 



30, 



(2a) Q 



and 



(2A. 0 2 - 1) log - 6\ Q 



q»a = Cj^l-y. 3 ^ (\ 2 ~1)£ (5. X log ^-3, J^-) 

^2 A.— 1 \ — 1 / 



CO S 0) 



Thus, the complete velocity potential consists of a part qf^ due to an axial velocity 
- RO cos 3, a part cp 2 due to a transverse velocity - Rf) sinj p\ and a part cp 3 due t>« a 
rotational motion of angular velocity 0 about OZi It follows then that 



(10) 



where 



n j.i v. „_„ ...j: <kj 



. . _ ... . ... * A n dp n dq n _ 



00 



d(Xk) 3P n dQ n d 

.TTTdirdT- ntI1+1) dH ^ 

- 2aftO sin p + p, ) as the boundary condition (reference 3) 



<P 3 ^lV Qn 1 ^) C0B 113 * ith Q =l B n J 

dp. 



and 



00 



cp 3 = n t 1 C n Pn 1 (M-) Qn 1 (X> cos u> with ^ C^ 1 



. dix dP dX n ^ J d|i^H a Jj - 



(2a) a n ((It A tM ^ 8 the boundary condition 
* dpi ' 

From the similarity of the left-hand side* of the boundary conditions assopiated with cp^, } 
it is clear that the treatment for cp 2 applies equally well to cp 3 . Thus, a method has 
been outlined for obtaining the velocity potential in a fluid due to a body of revolution in 
circular flight. 
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THE PRESSURE FORMULA 



The nest step is to calculate the forces acting on 
the moving "body owing to the pressure of the fluid medium. 
Since, the "body moves at a constant angle of attack and in 
a circular path with constant angular velocity, the veloc- 
ities V, t* are independent of the time so that the motion 
is a steady one. Therefore, the velocity potential cp 
does not contain the time t explicitly, i.e., = 0» 

and Lamb's formula (reference 1, p. 18) for the pressure 
when the coordinate axes are in motion becomes: 

(*g - -ID - -r (-g - 4?) - »° (*H - 'g) 

wnere co x , tOy, gj z are the component rotations referred to 

2 3 " 2 

the moving axes and q 3 = (urU) . + (v-V) + (w-'ff) where 
(u,v,w) = (- - |e. - ff). also cp"= <p x + cp 3+CP 3. " 

:If the radius vector to the point (x,y,z) is denoted hy 
r, the fluid-velocity vector "by - grad cp, and the rota- 
tion by U then equation (ll) may be written in vector 
notation as: 

f = - | (a-q) - (^£r^aa<?3) Ciaj 

where the symbols ( ) and [ 3 denote, respectively, the 
scalar and vector products of two vectors. The scalar 

triple product (£!■ [F- grad cp] ) = ( Cf5' ?] • "grad cp) = (t-grad cp) 
where t denotes that part of the velocity vector of any 
point attached to the moving axes due to the rotation 
alone. Equation (12), applied to the prohlem of this pa- 
per, then becomes: 

2. = - % (V 3 + V 3 ) - - ( grad cp . gFaT cp ) - + F* grad cp ) 

p 2 V <i ' ■ " _" i 

or, omitting ^ (V s 2 + Vy 3 ) ' since it is a constant and 

therefore has the same value over the surface of the "body, 
it follows that: 

p = \ ^XJ^'k-yS^ " \ (£^JP + VJ^-^ra^P+^J^J") (13) 
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In order to interpret equation (13) in spheroidal coordi- 
nates, it is necessary to obtain the presentation of the 
various vectors involved, in the Cartesian frame attached 
to the point .(p.. X, U>) • The .table of direction cosines 
for the ■ Cartesian .frame attached to any point is: 



y ■ 



X -1 



z 



: r- 

. e ■ * 
■X -I 



P 



V X' -p 



( X -1 \ 

~P ( 72 2~ ) COS (0 



CO 



-sin CO 



t 



\ ~P y 
cos to 



sin oj 



(14) 



With the use of this table of direction cosines the vector 
q , which for the present problem has the components 

3<P 3<p 3cp • • 

" 3x ~ x ' ~ 3y ~ y* " 3z W ^*'' :1 respect to the instantane- 
ous position of the moving axes, has in the Cartesian frame 
attached to any point \x, X, go the following components: 



1p 



■ 39 



*oo = 



39 
5iw 



"P 



r 



(15) 



" Voo 



J 



where 



B 3L = JL Adskl-Y i2 9JL _ _i_ /2L-i_\ 



>\*-p 5 



^ 3 -p? 



39 



and 



3iL = 
3s, 



00 



2a (X 2 -1) 



(1-P ) 



39 
'3co 
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i 



S CO, 



cos co, 7(0 = *- V Tr sin 0) 



Similarly, the vector which in the moving OXYZ frame 

has the components (0 , 0, Q) , has in the (p,, \, co) Car- 
tesian frame the following presentation: 



a 



0) 



cos CO 



sin co 



(16) 



and the radius vector r, the following components: 



l 



r co = 0 



(17) 



The vector IT has the components - Qy , 0 along 

the axes of the- moving frame and the following components 
relative to the rectangular (jx, X, co) axes: 



t„ = - 2a fiX CV 1 -^ cos CO 



» cos CO 



t„ = - 2a QXjj, sin CO 



(18) 



1-0 
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THE LONGITUDINAL- COMPONENT 07- THE CENTRIFUGAL 
FOECE OF THE APPARENT MASS 



If .„ 9 is ..the inclination to. the . X .axis of the out- 
ward-drawn normal to' the meridian curve at the point P, 
then the longitudinal force is 

B 2 2TT • ■ • . - 

F T =■- -/. f . p cos .9 hdtods . ■■ . ■ - 

. B i . ,o ..... ... 

"here ds is an element of length along a meridian curve; 
s L , s 3 the values of , s at the fore and aft ends of the 

body, and h the radius of the section through P perpen- 
dicular to the axis' of symmetry OX. 

putting -. ■ 

■ dh 

- cos 9 = t— ■ 
a s 

it follows that: 

2TT ... - 

*l = r S 1 pawah 3 ' 

a c 0 

The area of the section of radius h is A = TTh 3 and the 
average pressure at this section is defined to be 

2 IT 

P- = i~ / P d W 

It then follows that 

F T = / P'dA (19) 

c 

taken over the generating contour of the body of revolu- 
tion. 
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THE TRANSVERSE COMPONENT OP THE CENTRIEUG-AI 
FORGE" 03? THE APPARENT MASS 



The component of the pressure p in the plane XY. 

is p cos co and the resultant transverse force is there- 
fore ....... " - ■-_ , tr - . - 

S s 2TT 

j? T = - / / p cos co sin 6 hdco ds 

s x o - 



or, since gin 0 = 

ds 



= - J J hp cos co dcodx 



The resultant of the components of the pressure p .in the 
XY plane taken. round a section of radius h perpendicu- 
lar to the X axis, is given "by 

P" = h / p cos co dco 

o ' 

acting normally to the' meridian curve in the XY plane. 
Then 

2 T = ■» / P" dx (20) 

THE YAWING MOMENT 

The yawing moment about the axis of Z is given hy 

% 

M z = / P"l ds 
s i 

where I is the perpendicular distance from the origin 0 
to the line of action of P" (or the normal at the point 
of the meridian curve in the XY plane' at which P" is 
applied). The direction cosines at a point P (x,y), of 
the outward-drawn normal to the meridian curve in the XY 
plane are ; 

dy dx 
ds ds 
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It follows, by means of -elementary analytic geometry, that 

xdr + ydy ■ dr 

I = - — -r — = — r t— 

ds ds 

and hence • • • •" " 

M _ = - | / ' pn d (x 3 + y 3 ) = - i- / P" dr s (21) 
2 c * c 

taken over the -upper half of the meridian profile in the 
XY plane. It is interesting to note that for a sphero, 

since r 3 = constant, VL Z = 0. The differential dr 3 may 
bo looked upon as a measure of the deviation of the "body 
from a sphero. 
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APPENDIX 
SOLUTION POR AN ELLIPSOID 



The caso of an ellipsoid of revolution moving in a 
circular path has already boon worked out, but the methods 
used are very cumbersome and involve unnecessary algebraic 
manipulations. It seems therefore instructive of the pres- 
ont method, which makes liberal use of vector analysis, to 
rodovolop the solution for an ellipsoid. 

The equation of tho, meridian ellipse in Oartosian co- 
ordinates is: 

+ K = i 

A B 

and in ellipsoidal coordinates, simply 

* = \> 

where 

\ 0 = J, A 3 -B 2 =(2a) , Aq=(2aj , and 2a (\ -l) 2 = B 



The velocity potential <+> is the sum of the following three components (reference 1, 
pp. 132 and 133)*; 



u =■ - V cos p and v = - V Bin p 

The pressure formula (equation (13)) is expressed in vector form and is therefore 
independent of the particular coordinate system used* The components of the vectors 
^rad cp, V, and t are given by equations (16) to (18). For the simple case of an ellip- 
tical profile they are as follows: 





> (22) 



where 



COB UJ 



grad cp+V+t = < 



0 " 

' \ 

- '|v(l+M)/* A Op, U+e 2 K)| Bin 



U) 



cob ia 



V + t E=< 



where 



cos yj 



- (,v + AJlfl) Bin CU 



1- log 1+e - 1 



2e 1-e 1-e' 



I log Ira _ b 

_ T 2 1-e 1-e _ . . 

— li, . ^ - - M| 

2 b l- e 1-e 3 



3. lQg l±e__ 3 _ 
2e 1-e 



1-e' 



3 2-e a _ 1+e , e a 

log- 6 5 

2 e ^1-e 1-e 3 



= - H 



(23) 



124) 



2T.A.C.A. Technical- 3Tote Ho .. 554 



15 



Calculation of the longitudinal force Component 

According to equation (19) 

JFL = / P ' dA = 17 / P ' cL^i 2 = - 2rrB 2 Z 1 P'p dp, 

C C _1 

where 

The pressure p is easily calculated by forming the sca- 
lar products of equation (13) "by means of the vector com- 
ponents of equations (23) and (24) . It is useful to ob- 
serve that only the coefficients of cos 2 to, sin 2 dl need 

2TT 2TT 

"be considered since / cos (Odco = / sin codco = tt 

2TT 0 2TT 0 . 

while / cos wdoo = / sin oodfc) =0. It then follows 

o o 
after some algebraic reductions, that 

? L = | n A B 2 p Q M v = ~ m y. sin p 

• a a 

where m = — tt A B p, the mass of the fluid displaced by 

the spheroid. Noting that U = kj , the transverse iner- 
tia coefficient of the body: 

F L = - m k T |- sin P (25) 
Calculation of the Transverse Force Component Frp 



According to equation (20) 

i 

% = - / P" dx = - A / P" dp, 

c -1 

where 

2TT 

P" = h / p cos (Odd) 
o 

Again, observing that only the coefficients of coso) in 
the expression for the pressure p need be considered, it 
follows that: 
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Fm = - | tt A B 2 p QLti = bil|~cos p 
Noting that L .'= k_ , the axial inertia coefficient of the 

& » ■ , 

body: 

P T = m k a J- cos 0 (26) 



Calculation of .the ■ Yawing Moment M z . 

According to . equation. (21). 

M z '= -J-/ E>-» d(x 2 + y s ) •' ' 
. c 

NOW ■ . . ■ i ■ X 

x = 2a \ jj, and y = 2a(A. 0 - 1)2 (l - |_l, s ) 2 

so that s i 

M z = « (2a) / P" [X dp, 
- i 

Again, only .the . coefficient s of cos u> in . the. express ion 
for the pressure p need be considered and it follows 
that : 

2 



Prom the expressions, for L , M it is easy to- see that 

+ _JL_ = i 

2(L+1) (M+l) 

or .... - 

^ + (1+L) = (L+l)(M+'l) 
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It may be verified that 

■ 3 M+l e 3 V2 ■ S 1-6 5^1 / 3 
It follows that 

M z = - m u v (1+L) (2M - 1) 
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and since 

(1 + L) (2M - 1) = M - L = k T - k a , 

that 

M z = - Sg.- (k T - k a ) sin 2(3 (27) 
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